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Abstract 

The space-time conservation element solu- 
tion element (CE/SE) method [1] is em- 
ployed to numerically study the near- field 
of a typical under-expanded jet. For the 
computed case — a circular jet with Mach 
number Mj = 1.19 — the shock-cell struc- 
ture is in good agreement with experimen- 
tal results [2. 3]. The computed noise field 
is in general agreement with the experi- 
ment. although further work is needed to 
properly close the screech feedback loop. 

1 Introduction 

An under-expanded supersonic jet radiates mixing 
noise, broadband shock-associated noise, as well as 
screech tones under certain conditions. These com- 
plicated and technologically important physical phe- 
nomena have been the topic of many experimental 
and theoretical investigations, see Tam's review pa- 
per [4] for a comprehensive list of references. Sim- 
ply put. the mixing noise is directly associated with 
large-scale structures, or instability waves, in the jet 
shear layer, whereas the broadband shock-associated 
noise and screech tones are associated with the inter- 
action of these waves with the shock-cell structure 
in the jet core. The distinct screech tones arise due 
to a feedback loop. i.e.. part of the acoustic waves 
generated by the wavc/shock-ccll interaction propa- 
gate upstream and re-generate the instability waves 
at. or in the vicinity of, the nozzle lip. More de- 
tails can be found in the review papers [4-7] and the 
references therein. Screech tones are of particular 
interest not only because of general noise-reduction 
concerns, but also because of potentially destruc- 
tive structural interact ion (s) leading to sonic fatigue. 
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The feedback loop leading to screech tones is very 
sensitive to small changes in the system conditions 
and the understanding of the phenomena is to date 
mostly based on experimental observations, e.g. [2, 
3]. A modular simulation tool has been developed 
at McDonnell Douglas Aerospace, with components 
approximating the relevant physics, allowing design- 
ers to reasonably well estimate jet screech tones, see 
[7-9] for further details. Recently. Shcn and Tam [10] 
obtained good results in a direct numerical simula- 
tion of screech for a (low) supersonic Mach number 
axisymmetric jet. 

Jet noise is a challenging topic in computational 
aeroacoustics. in particular, near- field noise compu- 
tation in the presence of shock cells in the jet core. In 
this situation, the computational scheme is required 
on one hand to resolve the acoustic waves without 
introducing too much dispersion error and numerical 
dissipation, while on the other hand, it is required to 
capture shocks, or other nonlinear phenomena, near 
or inside the jot correctly. In addition, non-reflecting 
boundary conditions must be implemented, which is 
more difficult to accomplish in the near field than in 
the far field. 

The recent Space-Time Conservation Element 
and Solution Element Method 1 [11]. or the CE/SE 
method in short, is a scheme that meets the above re- 
quirements. The CE/SE scheme possesses attractive 
properties for aeroacoustics computations in that: 
(i ) it possesses low dispersion and dissipation errors; 
( v ) its ‘built-in’ shock-capturing nature makes the 
computation of shock-cell structures simple and ac- 
curate: (Hi) the non-rcflccting boundary conditions 
arc simple and effective and can be applied in the 
near field of the jet without introducing excessive 
errors: and (tv) the scheme delivers the design ac- 
curacy for the vort icily, which plays an important 
role in the noise generating mechanism. A detailed 
description of the CE/SE method can be found in 
the reports of Chang ct al [12,13]. As demonstrated 
in our previous papers, the CE/SE scheme is well 
suited for computing waves on compressible shear 
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flows [14] as well as vorticitv /shock interactions [14, 
15], both being corner stones of the jet-noise phe- 
nomena. 

In this paper, the near-field noise of a shock- 
containing axisyinmetric supersonic jet. is investi- 
gated numerically by using a CE/SE Euler solver. 
The paper is arranged as follows: The axisyinmet- 
ric CE/SE scheme is briefly discussed in Section 2. 
Section 3 illustrates the initial and boundary con- 
ditions for the jet-noise problem, in particular, the 
novel non-reflecting boundary conditions which are 
based on flux balance. The numerical results are pre- 
sented and compared to experimental findings [ 2 , 3 ] 
in Section 4. Conclusions are given in Section 3 . 

2 Axisvmmetric CE/SE Euler Method 

In general, the CE/SE method systematically solves 
a set of integral equations derived directly from 
the physical conservation laws. Because of its in- 
tegral formulation and since, as will be seen below, 
space and time are treated ’on the same footing’, 
the scheme naturally captures shocks and other dis- 
continuities in the flow. Both dependent variables 
and their derivatives are solved for simultaneously 
and. consequently, the flow vorticitv can be ob- 
tained without reduction in accuracy. Non- reflect ive 
boundary conditions are also easily implemented be- 
cause of the flux-conservation formulation. 

2,1 Conservation Form of the Unsteady 
Axisymmetric Euler Equations 

Consider a dimensionless conservation form of the 
unsteady axisymmetric Euler equations of a perfect 
gas. Let p, u. v. p, and 7 be the density, stream- 
wise velocity component, radial velocity component, 
static pressure, and constant specific heat ratio, re- 
spectively. The axisymmetric Euler equations then 
can be written in the following vector form: 

Ut-+-F x + G,j = Q. ( 1 ) 

where x. y ^ 0 . and t are the streamwise and radial 
coordinates and time, respect ivcly. and the conser- 
vative flow variable vector U and the flux vectors in 
the streamwise and radial directions. F and G, are 
given by: 


c: 

11 

11 

( Gl \ 

c, 

Isi 

with 



£ 

II 

0 

5 

II 

f 3 = pc- 


U\ = p/( 7 - 1 ) - p(u 2 - V 2 )/ 2 


Fi = U-2 , 


F, = (-. - 1 )Ui + [(3 - -)Ui - (- - 1 )U 2 ] /2 U v . 
F 3 = U2V3/V1. 


Ft = 'UiUi/lh - (* - 1 ) U-2 \Ui + I' 2 ] / 2 V 2 . 

Cl =13. G 2 = t/ 2 l/ 3 /l/i, 

Ch = (7 - DC.1 + [(3 - - )l"i - - 1)C|] /2Ui. 

Cl = -IWi/ih - (- - 1 )U 3 [U 2 + V 2 ] /2 Uf, 


and 


with 


Q = 



\ (h 1 

\qJ 


Q 1 = -v 3 /y. q 2 = -v 2 u 3 /u iy . 


Ch = -Ui/U L y, ch = -Ci/i/. 

The axisymmetric Euler equations can be though 
of in this formulation as a variation of their two- 
dimensional counter parts with 'source’ terms on 
their right-hand sides. 

By considering (x. y, t) as coordinates of a three- 
dimensional Euclidean space £3 and using Gauss* 
divergence theorem, it follows that Eq. ( 1 ) is equiv- 
alent to the following integral conservation law: 



H m - dS = 



m = l,2.3,4. ( 2 ) 


where S(V\) denotes the surface around a volume \ * 
in £3 and H m = (F m .G m .U m ). 

2.2 CE/SE Structure 

In the CE/SE scheme, the flux conservation relation 
in space-time is the only mechanism that transfers 
information between node points. The conservation 
element CE . or computational cell, is the finite vol- 
ume to which the integral flux condition ( 2 ) is to be 
applied. Discontinuities are allowed to occur in the 
interior of a conservation element as long as the inte- 
grals on the right-hand side of (2) exist. A solution 
element SE associated with a grid node is here a sot 
of three interface planes in £3 that passes through 
this node. The solution, i.e. 17, U x , and U y arc 
calculated at this node. Within a given solution cl- 
ement SE(j.k, n). where (j, k. >?) is the node index, 
the flow variables arc not only considered continuous 
but arc also approximated by linear Taylor expan- 
sions: 


U'(x,y.PJ. k,n) = U'J k + (U^(x - x**)- 
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(U u yi,(v-yj.,) + (U t )l k (t-n. (3) 

F*(x,y. t;j. k.n) = F] k . + (F x )l k (x - Xj , k W 

(F <j)j.k(y ~ Vj.k) ~ (F tYj.kXt - t n )' ( 4 ) 

G*(x.y. t-j. k, n) = G'j k + (G x )J k (x - xj. k )~ 

(G y yi k (y-y M -) + {G t )” k (t-t n ). (3) 

where the partial derivatives of F and G can be 
related to the corresponding one of U by using the 
chain rule and Uf can be obtained from (1). 

The discrete approximation of (2) is then 

H* m • (IS = 1 '(CE) £ $j,k(Qm)j.k- 

{j,kji)<=S(CE) 

( 6 ) 

for in = 1.2,3. where H* m = (F^.G^.L^,) and 
E(j.k.n)esiCE) = !• The right-hand side of (6). 
in general* is the volume of the CE times a weighted 
average of the ‘source’ term evaluated at the nodes 
on S{CE). Each S(CE) is made up by surface seg- 
ments belonging to two neighboring SE's. All the 
unknowns are solved for based on these relations. 
No extrapolations (interpolations) across a stencil 
of cells arc needed or allowed. 

In order that the number of equations derived 
from the flux conservation law above matches the 
number of unknowns (here 12 scalar unknowns), the 
grid needs to be carefully designed [11, 13]. The 
mesh is staggered in both time and space. In a spa- 
tial plane in E3, the grid nodes, see Fig. 1, are 
grouped as two staggered sets Hi (open circles) and 
fb (filled circles). At a given time step the unknowns 
are evaluated only on one of the grid sets, fii or <>2- 
and at the next, it will be evaluated on the other set. 
i.c...the spatial sets alternate as time is stepped for- 
ward. As can be seen in this figure, each 'interior’ 
node point then has three nearest neighbors at the 
previous time step. Thus, there arc three CE' s as- 
sociated with each node point in this arrangement 
and. therefore, there arc the same number of rela- 
tions as there are unknowns. Fig. 2 illustrates the 
three conservation elements associated with a node 
point and its three nearest neighbors at the previous 
time step. Because of the staggered spatial grids, the 
method can be thought of as a 'two-step' method in 
time. 

2.3 Treatment of the Source Term 

In solving (6), only the 3’j k corresponding to the 
new (half) time level is taken to be nonzero (and 
equal to unity). This strategy has successfully been 
applied by Yu ct al [16] to problems with very stiff 
source terms. 


i 


S(CE) 


As the source term Q = Q(U) itself is a func- 
tion of the unknown U at the new (half) time level, 
a local iterative procedure is needed to determine 
U. For our particular choice of 3"*, the discretized 
integral equation (6) reduces to the form 

U-Q(U)^ = U„. (7) 

where U h is the local homogeneous solution (Q = 0 
locally). Note that U h only depends on the solution 
at. the previous (half) time step, i.e. U h is obtained 
using explicit formulas. A Newton iterative proce- 
dure to determine U is then 

t/( ,+1) = u {i) - 

where i is the iteration number and 
*(U) = U-Q{U)3L. 


Normally, U at. the previous (full) time step is a good 
initial guess U^ 0) and the procedure takes about 2-3 
iterations to convergence. The Jacobian matrix is 
given by 
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At 7-H/! + 3f7i , -17.1- 
2 y [ 2 Uj + Ui J 


The inverse of the Jacobian, i.e., (^§) _A can easily 
be derived analytically for this particular case, thus, 
leading to savings in CPU time. 


3 The Jet-Noise Problem 

Consider a circular jet as sketched in Fig. 3. The 
flow inside the sonic jet nozzle is a choked flow with 
exit Mach number M, = 1, while outside the noz- 
zle in the jet. the Mach number is Mj — 1.19. and 
the ambient flow around the jet is stationary. These 
conditions correspond to the experimental setup of 
Panda [2. 3]. In his experiments, it was shown that 
for the above Mj = 1.19. the jet. noise field is in 
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an axisymmctric mode. Hence, rhc axisymmctric 
CE/SE code described above is appropriate for com- 
putation at these conditions. 

In this investigation, our attention is focussed on 
the near field of the jet nozzle since this is the source 
region for acroacoustic noise. The diameter D of 
the jet nozzle is chosen as the length scale. The 
computational domain is a circular cylinder of Y2D 
axial length and oD radius. The grid currently used 
is a uniform one and the numbers of nodes in the x 
and y directions are 301 and 401, respectively, since 
good grid resolution is needed in the jet shear layer. 
The computational domain is about the same size 
as the experimental one in the radial direction but 
it is about twice as long in the axial direction. The 
computational scheme is of the a-e type [12,13] with 
o = 0 and e = 0.5. 

3.1 General Initial Conditions 

As shown in Fig. 3. there are four regions in the ex- 
periment. namely, the plenum, the convergent noz- 
zle with choked flow, the jet core area, and the sur- 
rounding ambient flow area. Only the latter two 
regions (III and IV in the figure) are included in the 
computational domain. A conceivable initial condi- 
tion is to assume that the pressure is uniform in the 
entire computational domain and equal to the am- 
bient pressure, except at the jet nozzle exit where a 
higher pressure p, is imposed since the jet is under- 
expanded. We will also make the reasonable assump- 
tion that the temperature in the plenum equals the 
ambient one. 

Let the density, axial velocity, and temperature in 
the jet core, be the scales for the dependent vari- 
ables. Hence, the initial conditions in the jet core 
are given by 

Pj = L T J = 1 • u j = !• V J = 0- Pj = zjjl- 

where T denotes the temperature. The initial con- 
ditions in the ambient region are given by 

Pa = \/T n . T„ = 1 + £(7 - l)iVj, 

U a - 0. V a = 0. P„ = ZJJ2- 

where the result for the temperature follows from 
the assumptions of constant total enthalpy in the jet 
flow and that the plenum and ambient temperatures 
are equal. 

By using the assumption of constant total en- 
thalpy and the condition for iscntropic flow, it fol- 
lows that the pressure p c at the nozzle exit can be 


related to the pressure p; in the jet core, thus, yield- 
ing 

P( _ + 

Pj [l + 4(7-l)J/2_ 

3.2 Initial Shear Layer Assumptions 

The ambient and jet core regions are connected by 
a continuously changing shear layer. In the present 
numerical computations, two distinct types of ini- 
tial shear layers were assumed and it was found, as 
will be seen below, that the results were somewhat 
sensitive to this choice. 

3 .2 A Parallel Shear Layer The shear layer is 
assumed to be of constant thickness and is taken 
to be a hyperbolic- tangent one. 

u(y) = |[1 - tank 20(</ - 5)]. 

With the Prandtl number taken as unity for simplic- 
ity, the temperature in the shear layer is then given 
by the Buseman-Crocco’s relation 

T(y) = T„ + (1 - T a )u - i(- - - u), 

which in this case reduces to 

ny) = l + kh-l)Mj{l-u a ). 

The density profile in the shear layer is simply given 
by p(y) = 1 /T{y). 

3.2.2 Artificially Expanded Shear Layer The 
empirical quasi-parallel shear-layer profile described 
by Morris and Tam [17] for a fully expanded jet is 
imposed as the initial condition. This shear layer has 
a constant pressure and zero radial velocity, but the 
axial velocity depends weakly on the axial coordi- 
nate through slow changes in the mixing layer width 
b, the radius of the potential core h, and. ultimately, 
the core velocity u c . 

u = u c . y ^ h: u = u c exp[- In 2 ( - - - ) 2 ]. y ^ /;; 

where 

1.35 

u c = 1, x ^ x c : u c = 1 - exp(- — — — ), x ^ x c ; 

1 — x/x c 

b = 1 “^ 0 _ 8 - r (i _ .273 Mr). X c = 4.2 - 1.1 Mj. 

x c is the length of the potential core. The radius 
of the potential core, h . is determined by consider- 
ing the conservation of axial momentum. After u is 
determined, the density p can be obtained as above. 
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4 Numerical Results 


3.3 Boundary Conditions 

At the inlet boundary, the dependent variables and 
their spatial derivatives are specified to be those of 
the ambient flow, except at the nozzle exit, 0 ^ y ^ 
where the elevated pressure p — p c is imposed. 
i.e. the jet is under-expanded, as in the physical 
experiment. 

An artificial forcing of the shear layer is also in- 
troduced at the lip of the nozzle through a small 
perturbation in the radial velocity component. 

r'(x = 0) = <>sin(27rS<)cxp[— (y — ±) 2 ], 

in order to organize the flow somewhat and to pro- 
vide a reference signal for conditional sampling of 
the computed results. The amplitude 6 = 0.001 is. in 
general, used throughout this study and the Strouhal 
number, or nondimcnsional frequency. St = 0.59 
corresponds to the experimentally observed screech 
tone. 

At y — 0. the appropriate conditions are applied 
which can be deduced from the limiting form of 
Eq. (1) as y — ► 0. At the top and outflow bound- 
aries. the Type I and Type II CE/SE non-reflecting 
boundary conditions as described in the next sub- 
section arc imposed respectively. 

3.4 Non-Reflecting Boundary Conditions 

In the CE/SE scheme, non-rcflectivc boundary con- 
ditions arc constructed so as to allow fluxes from the 
interior domain to a boundary CE smoothly exit to 
the exterior of the domain. There are various vari- 
ants of the non-reflecting boundary condition and in 
general they have proven to be well suited for acroa- 
eoustic problems [18-23]. The following are the ones 
employed in this paper. 

For a grid node (j.k,n) lying at the outer radius 
of the domain the non-reflcctive boundary condition 
(type I) requires that 

(U. v )” k = (U v %> = 0 . 

while U n j k is kept fixed at the initially given steady 
boundary value. At the downstream boundary, 
where there arc substantial gradients in the radial di- 
rection. the non-rcflectivc boundary condition (type 
II) requires that 


(U x ) » k = G. 

while U n f k and (U ,j)j k are now defined by simple 
extrapolation from the interior, i.e.. 


U 


n 

j.k 


= u 


n- 1/2 

jb*' 


(U 


yfj,k 


= (U 


\n- 1/2 
u)j\k' 


In this section, the numerical results for the axisym- 
mctric under-expanded circular jet are presented 
and compared to experimental results [2. 3]. At 
the low supersonic jet Mach number Mj = 1.19 
and the exit pressure to ambient pressure ratio 
p ( /pa = 1.2645, the overall motion in the exper- 
iment [2] is in an axisymmctric mode and with a 
screech tone at the Strouhal number St — 0.59. 
For imperfectly expanded jets such as this one. a 
quasi-periodic shock-cell structure is formed in the 
jet plume. The strcamwisc-growing instability waves 
naturally occurring in the jet shear layer interacts 
with this shock-cell structure and generate broad- 
band shock-associated noise and under certain con- 
ditions screech tones through a feedback loop. 

It is important to realize that even if no harmonic 
forcing is intentionally imposed on or introduced in 
the numerical simulation, the jet shear layer is in ac- 
tuality continuously stimulated at a very low level as 
a result of truncation, round-off. and discretization 
errors (all of which can be characterized as numeri- 
cal noise — and. in a sense, are analog to environmen- 
tal background noise in experiments). Among these 
(random) perturbations, those whose frequency falls 
within the unstable range for the jet shear layer will 
draw energy from the mean flow and grow in the 
streamwise direction, and the other will eventually 
decay out. For the initially parallel shear layer above 
at the conditions under consideration, the unstable 
range is 0 < St < 0.6963. These growing perturba- 
tions interact with the shock-cell structure of the jet 
plume and thereby generate acoustic waves. 

As has already been pointed out, a small (ampli- 
tude 10“ 3 ) harmonic forcing at the observed screech 
Strouhal number is imposed on the shear layer at 
the nozzle lip. where its receptivity is considered the 
highest, in our computations. This will slightly or- 
ganize the flow and thereby allow conditional sam- 
pling (or phase averaging) techniques to be used to 
extract wave generation and propagation properties 
relevant to the screech phenomenon from the nu- 
merical computation. In an ideal simulation, the 
harmonic forcing could be turned off once the self- 
sustained screech tone has developed through the 
feedback loop, or might not even be necessary. How- 
ever. in our simulation, the nozzle lip is of zero thick- 
ness and therefore' we do not expect the feedback 
loop to be 'properly closed' at present. This expec- 
tation is borne out by the computations and a finite 
nozzle lip will be used in planned future work. 
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4.1 Jet Shock-Cell Structure 

Experimental results for supersonic jets are often 
documented in terms of Schlicrcn pictures. It is 
straightforward to construct Schlicrcn graphs from 
the numerical results without loss of accuracy in 
view of the nature of the CE/SE method. Fig. 4 
shows a comparison of numerical Schlicrcn (density 
gradient modulus) contours from the current com- 
putations for the initially parallel shear layer case 
and the experimental Schlicrcn photograph in [2]. 
The numerical result is a snapshot of the instan- 
taneous density gradient after about 6 periods of 
computation from actual startup of the simulation, 
where a the period is defined in terms of the imposed 

Strouhal number, i.e. 1 /St = 1.6949 Due to the 

nature of the Euler solver, there is practically no dis- 
sipation in the solution, thus strong vortex roll-up 
occurs in the shear layers downstream of the shock 
cells in contrast to the blurred shock cells observed 
in the experiment. However, the shock cells that 
have developed so far in the simulation agree quite 
well with the experimental ones. 

4.2 Near-Field Radiating Waves 

As have already been alluded to, this highly nonlin- 
ear system with discontinuities (shock-cell structure) 
exhibits wave motion in a rather random fashion. In 
order that the waves associated with the particu- 
lar frequency of St = 0.59 can be clearly displayed 
and examined, a conditional sampling technique is 
applied to the computed data. Using the harmonic 
forcing as a phase reference, the computed results 
can be averaged over several periods for any partic- 
ular phase angle selected. This technique is some- 
times referred to as phase averaging. In the present 
work, this procedure was carried out for 24 phase an- 
gles using 50 cycles to construct the averages after 
the start-up transients propagate out of the compu- 
tational domain (12-15 cycles). It should be noted 
that the conditional sampling (or phase averaging) 
for efficiency is implemented as an integral part of 
the time marching numerical procedure. 

Fig. 5 demonstrates a typical result (one of the 24 
phase averages) for the initially parallel shear layer 
case. It is a combination of isobar contours and nu- 
merical Schlicrcn (moduli of density gradients) con- 
tours. In the processing of the numerical data for 
this figure, the very high level isobar contours, corre- 
sponding to hydrodynamic waves around the jet core 
area are ‘cut off and the ‘colors’ arc appropriately 
adjusted so that the acoustic waves are clearly dis- 
played. As seen in Fig. 5, these acoustic waves prop- 
agate in both upstream and downstream directions 
separated by a radial line. The first harmonic — 


waves with half the period or twice the frequency — 
is expected to be dominant along this dividing line. 
The acoustic wavelength is about 1.6 D. each shock 
cell (a pair of shock diamonds ) spans about 0.8D 
in the axial direction. These results agree well with 
experimental results (c.g. [2]) and theoretical pre- 
dictions (e.g. Seiner [5]. Tam [4]). 

The remaining phase averages are not included 
here for brevity, but all 24 can be seen as an anima- 
tion movie at the web address 1 

www.grc.nasa.gov/WWW/microbus/ccsc/scrccch html 

In the animation, it is observed that the first shock 
cell does not appear move or deform at all, i.e. the 
shear-layer instability wave is too weak at this lo- 
cation to significantly affect the shock cell. How- 
ever. it interacts strongly with the shock cells once 
it has achieved a sufficient amplitude through its 
strcamwisc growth, which is evident in the anima- 
tion by the motions an deformation of the third and 
fourth shock cells. The deformation and movement 
of the shock cells arc also seen in the experiments 
[2. 3]. This is, as pointed out by Tam [4]. the 
mechanism by which the acoustic waves involved in 
broadband shock-associated noise and screech tones 
arc generated and this is also evident in the ani- 
mation. Furthermore, the fifth shock coll and those 
further downstream in the animation are so severely 
deformed that they appear quite different from tiie 
first few shock cells. 

Fig. 6 is a side-by-side comparison of the com- 
puted acoustic field for the case with the artificially 
spread shear layer and the experimental one [2] for 
four phase angles. Note that for this comparison, 
only half of the computational domain (0 ^ x ^ 6D) 
is shown. Both the experimental and computational 
data in this figure arc phase averaged results. The 
increment in the phase angle between the different 
panels is the same for both type of results. How- 
ever. it is not possible to assign an absolute phase 
angle to both cases since the phase difference be- 
tween the ‘reference’ signals used to obtain the con- 
ditional samples in the two cases is not known. Con- 
sequently. the panels arc arranged in such a way 
that their phases in a period match the best to each 
other. It can be observed that the upstream prop- 
agating acoustic waves match well with the experi- 
mental ones, while the numerical downstream prop- 
agation waves appear to be too strong. 

Fig. 7 shows the phase averaged pressure fluctua- 

~A frames version is also available at 
http://www.grc.nasa.gov/WWW/microbus, which gives 
access to further CE/SE examples and information. 
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tions at the nozzle lip at various phase angles during 
a period. These waves arc obviously non-sinusoidal. 
they contain other frequencies even though the forc- 
ing at the nozzle lip is a low-level pure sine wave. 
Thus, there is feedback to the nozzle lip from the 
strong shock-cell/instability- wave interactions oc- 
curring further downstream. 

A simple test was conducted by turning olf the 
artificial perturbation of the shear layer at the jet 
nozzle after about 50 periods and running an ad- 
ditional ten periods or more. The overall level of 
the unsteady motion remained about the same, the 
motion became less organized, and was judged dom- 
inated by the natural instability of the shear layer. 
This implies that the feedback loop is not properly 
closed at the lip in the present simulation. 

4.3 Time- Averaged Flow 

As currently implemented, the conditional sampling- 
procedure provides averaged results corresponding 
to 24 phase angles. The time average results are 
obtained by averaging over these 24 phase angles. 

Fig. 8 -depicts time-averaged results for the case 
of an initially parallel shear layer. Part (a) shows 
the axial velocity u profiles as well as contours. 
Spreading of the shear layer due to mixing is evident. 
Part (b) shows the pressure and density distributions 
along the center line (y = 0). The pressure and den- 
sity distributions arc in qualitative agreement with 
experimental results, included here for convenience 
in Parts (c) [5] and (d) [2]. 

Fig. 9 shows the corresponding time-averaged nu- 
merical results for the case of the initially artifi- 
cially spread shear layer in Parts (a) and (b). The 
mean shear layer here spreads more than in the 
previous case, but, somewhat surprisingly, it ac- 
tually spreads slower than prescribed by its initial 
condition. Please recall, however, that the initial 
condition used corresponds to a perfectly expanded 
jet and. obviously, the under-expanded situation of 
the actual computation then must lead to a re- 
adjustment of the mean shear-layer flow. Also note 
that the mean pressure and density variations along 
the center axis attenuate faster than in the previous 
case. 

Fig. 10 compares time-averaged Schlicrcn graphs 
corresponding to the two initial conditions for the 
jet shear layer. Obviously, the shock cell structure 
in the second case decay faster, but look closer to the 
experimental one in Fig. 4 than the first one docs. 
The sonic line within the mean shear layer slants 
slightly more upwards in the second case, allowing 
the shock cell structure to decay faster, while in the 
first case, the sonic line is more parallel to the center 


line thus causing the annular shear layer to act as 
a better waveguide and. consequently, the shock-cell 
structure survives longer. 

4.4 Power Spectra 

The numerical time history is recored for a select 
number of locations in the flow field and later post- 
processed to obtain spectral information using Fast 
Fourier Transform (FFT) techniques. Estimates of 
the power spectral density (PSD) for pressure arc 
shown in Fig. 11. where Part (a) is for the case 
with the initially parallel shear layer and Part (b) is 
for the initially artificially spread shear layer. The 
Figure shows results for three locations at approx- 
imately 5 D distance from the origin (taken to be 
located at the center axis at the jet exit) for each 
case. The solid and long-dashed curves are for a po- 
lar angle of 61° and 29°. respectively, in both cases. 
The short-dashed curves arc for a polar angle of 11° 
in Part (a) and 5° in Part (b). These power spec- 
tral density estimates arc typical ones for the flow 
field. Because of the comparatively short time pe- 
riod (about 50 periods) available for the spectral 
analysis, a trade-off has to be made between the 
spectral resolution and the reduction in variance of 
the estimate. Even though somewhat crude, the 
results shown represents the best balance between 
these two conflicting requirements. 

The short-dashed curves in Fig. 11 represent 
power spectral densities at or near the shear lay- 
ers. The forcing Strouhal number {St = 0.59) is dis- 
cernible in these curves, but it is clear that at these 
locations the unsteady motion is dominated by the 
natural instability waves of the jet shear layer. The 
overall level is the highest for these curves since the 
strong hydrodynamic waves arc involved. As the 
’observation' point moves away from the vicinity of 
the shear layer, the overall level will be reduced as 
the unsteady motion becomes dominated by the gen- 
erated acoustic waves. 

The solid curves in Fig. 11 correspond to a lo- 
cation in the upper left corner of the domain, 
where acoustic waves dominate the unsteady mo- 
tion. There arc here substantial differences for the 
two cases. In Part (a) the curve shows that the 
mixing noise is dominating the unsteady motion, 
whereas in Part (b) the curve shows a clear peak at 
the forcing frequency and significant contributions 
to the acoustic field by broadband shock-associated 
noise. The mixing- noise level is also somewhat 
higher in Part (b). It is clear that due to the dif- 
ference in the shock cell structure (Fig. 10). par- 
ticularly we believe due to the more rapid decay 
of the shock-cell structures further downstream, a 
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much stronger noise field is generated in the second 
case. 

5 Concluding Remarks 

In this paper, the recent CE/SE Euler scheme is ap- 
plied to noise computation for a supersonic shock- 
containing axisymmctric jet. The advantages of the 
CE/SE scheme were confirmed in the present sim- 
ulation. The implementation is ‘effortless' in that 
no special treatment and parameter selections are 
needed. Many aspects of the computed results are 
also in good agreement with experimental findings 
[2, 3]. However, from the literature [4] and our expe- 
rience herein, it is clear that both the annular shear 
layer spreading near the jet nozzle and the details of 
shock-cell structure affect the overall quality of the 
noise prediction. It is concluded that the sensitivity 
of our results to the initial shear-layer assumptions 
is inherent in using an Euler scheme. We plan to ap- 
ply a CE/SE Xavier-Stokes solver, desirably with a 
large-eddy simulation capability in our further inves- 
tigation. There is reason to believe that this should 
not only provide a correctly spreading jet shear layer 
but should also reduce the amplitude' of the down- 
stream propagating acoustic waves. Finally, a self- 
sustained oscillation, i.e. the screech tone, is yet to 
be achieved in our simulation. The addition of a fi- 
nite lip of the nozzle coupled with the viscous scheme 
should rectify this situation. The results will be re- 
ported in the future. 

The authors wish to thank Dr. J. Panda for fruit- 
ful discussions about the jet screech phenomenon. 
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Figure 1: The spatial node locations of the grid sets 
fii (hollow circles) and (solid circles), showing 
the hexagon structure of CE (1} = CEj 1 ^ — CE.^ — 
CE^ with STi ,; shaded hexagon corresponding to 
Fig. 2(a) and x-axis being the axis of symmetry. 
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Figure 2: (a) Conservation clement CEj^ and (b) 
solution element SE{^ 


12D 



Figure 3: Sketch of the jet noise problem: only Re- 
gions III and IX are included in the computational 
domain. 
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Figure 4: Comparison of numerical Schlicrcn graph (instantaneous snapshot, left panel) and the experimental 
time averaged Schlicrcn graph (right panel); showing excellent agreement of the first few shock cells. 


radiating 

waves i «n 



12D 


MJ=1.19 axisym. jet (underexpanded) 


Figure o: Typical phase averaged isobars and numerical Schlicrcn pictures (parallel initial shear layer) 
showing radiating waves and the shock-cell structure. 
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Figure 6: Comparison of upstream propagating wave patterns with experimental screech waves for phase 
phase angles (artificially spread initial shear layer). 
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time averaged u contours and profiles (simple shear layer) 




time averaged p and rho at axis 
(simple shear layer) 




Figure 8: Time-averaged u profiles and contours; time averaged p and p along x-axis; experimental and 
computed result by others (parallel initial shear layer). 
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time-averaged u, showing shear layer spreading 



Figure 9: Time-averaged u profiles and contours; time averaged and p along x-axis (artificially spread 
initial shear layer). 
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Figure 10: Time-averaged numerical Schlieren pictures, showing shock cell structures: (a) parallel initial 
shear layer: (b) artificially spread initial shear layer. 
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Figure 11: Power spectra! density (PSD) for the pressure versus Strouhal number: (a) parallel initial shear 
layer, (b) artificially spread initial shear layer; showing mixing noise and broadband shock-associated noise. 


16 

American Institute of Aeronautics and Astronautics 




REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden tor this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services. Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway. Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget. Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1. AGENCY USE ONLY {Leave blank) 


TITLE AND SUBTITLE 


REPORT DATE 

December 1999 


REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


Noise Computation of a Shock-Containing Supersonic Axisymmetric Jet 
by the CE/SE Method 


6. AUTHOR(S) 

Ching. Y. Loh. Lennart S. Hultgren, Sin-Chung Chang, and 
Philip C.E. Jorgenson 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
John H. Glenn Research Center at Lewis Field 
Cleveland. Ohio 44135-3191 


WU-522-8 1-1 1-00 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E- 12029 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TM— 1999-209658 
AIAA-2 000-04 75 


11. SUPPLEMENTARY NOTES 


Prepared for the 38th Aerospace Sciences Meeting and Exhibit sponsored by the American Institute of Aeronautics and 
Astronautics, Reno, Nevada, January 10-13, 2000. Ching Y. Loh, Taitech Inc., Cleveland, Ohio 44135: Lennart S. 
Hultgren, Sin-Chung Chang, and Philip C.E. Jorgenson, NASA Glenn Research Center. Responsible person, Philip C.E. 
Jorgenson, organization code 5880, (216) 433-5386. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


Unclassified - Unlimited 
Subject Category: 71 


Distribution: Nonstandard 


This publication is available from the NASA Center for AeroSpace information. (301) 621-0390. 


13. ABSTRACT (Maximum 200 words) 

The space-time conservation element solution element (CE/SE) method is employed to numerically study the near-field 
of a typical under-expanded jet. For the computed case — a circular jet with Mach number Mj = ] : 19 — the shock-cell 
structure is in good agreement with experimental results. The computed noise field is in general agreement with the 
experiment, although further work is needed to properly close the screech feedback loop. 


14. SUBJECT TERMS 


Space-time conservation element and solution element method, Aeroacoustics, Space- 
time flux conservation 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


15. NUMBER OF PAGES 

21 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 
















